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We present a summary of the Fermi Pulsar Search Consortium (PSC), an international collaboration of radio 
astronomers and members of the Large Area Telescope (LAT) collaboration, whose goal is to organize radio 
follow-up observations of Fermi pulsars and pulsar candidates among the LAT 7-ray source population. The 
PSC includes pulsar observers with expertise using the world's largest radio telescopes that together cover the 
full sky. We have performed very deep observations of all 35 pulsars discovered in blind frequency searches of 
the LAT data, resulting in the discovery of radio pulsations from four of them. We have also searched over 
300 LAT 7-ray sources that do not have strong associations with known 7-ray emitting source classes and have 
pulsar-like spectra and variability characteristics. These searches have led to the discovery of a total of 43 
new radio millisecond pulsars (MSPs) and four normal pulsars. These discoveries greatly increase the known 
population of MSPs in the Galactic disk, more than double the known population of so-called 'black widow' 
pulsars, and contain many promising candidates for inclusion in pulsar timing arrays. 
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1. The Fermi Pulsar Search Consortium 

The Fermi Pulsar Search Consortium (PSC) was es- 
tablished as a collaborative effort between the Fermi 
LAT team and an international group of radio pulsar 
observers with expertise using radio telescopes around 
the world. The two goals of the PSC are: (1) to ob- 
tain radio observations of 7-ray pulsars discovered in 
blind searches of LAT data, (2) to search for radio 
pulsars in the directions of unassociated LAT sources 
that have pulsar-like characteristics. This collabora- 
tion has proven extremely successful, as we describe 
below. Members of the PSC have used the following 
radio telescopes for these observations: Green Bank 
Telescope (GBT), Parkes, Effelsberg, Nangay Radio 
Telescope (NRT), Arecibo Telescope, the Lovell Tele- 
scope at Jodrell Bank, and the Giant Metrewave Ra- 
dio Telescope (GMRT). 

2. Radio Searches of Blind Search 
Pulsars 

One of the great successes of the Fermi LAT has 
been the discovery of 35 pulsars in blind searches of 
the 7-ray data alone [3 [2H GH E3] • Although these 
discoveries were made using the 7-ray data alone, be- 
fore these new pulsars can be dubbed "radio quiet", 
they must first be searched deeply for radio pulsa- 
tions. In addition, any radio detections provide im- 
portant new information on each pulsar including a 
distance determination from the dispersion measure 
(DM), which is critical for converting measured 7-ray 
fluxes into luminosities and thus determining the 7- 
ray efficiency of the pulsar. A radio pulse detection 
also yields information on the location of the emission 
region from the radio to 7-ray offset [3l] . Radio po- 
larization studies also provide potentially important 
information on the pulsar geometry (e.g. the angles 
between the rotation and magnetic axes [a] and be- 
tween the rotation axis and the line-of-sight [£]) [25] . 
Radio limits and detections also provide critical input 
for population models, where the blind search pulsars 
give a much more complete sample of the young pulsar 
population, owing to the broad 7-ray beams [351 I36j . 

The PSC has made extremely deep radio observa- 
tions of all of the LAT-discovered blind search pul- 
sars and detected radio pulsations from four of them 
[31 [HI [21] ■ Now, with 4 detections and 31 upper lim- 
its, it is appropriate to take a look at our progress so 
far. In Figure[T] we show the flux limits and detections 
(scaled to 1400 MHz using a typical pulsar spectral in- 
dex of —1.6) from our observations. It is clear that we 
are achieving flux limits competitive with the faintest 
known radio pulsars (all well under 0.1 mJy) and two 
of our detections (PSR J1907+0602 and J0106+4855) 
are very faint indeed! We also note that a possible 
(non-PSC) radio detection of PSR J1732-3131 has 
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Figure 1: Radio flux density upper limits (red arrows) 
and detections (green circles) for the 35 LAT blind search 
pulsars, scaled to 1400 MHz using a spectral index of 
— 1.6. The 1400 MHz flux densities of young pulsars in 
the ATNF catalog are plotted as black dots, for 
comparison. 

been reported at 34.5 MHz [TS] and so additional 
searches at low frequency are certainly warranted. 

The vast majority of these pulsars are so faint in the 
radio that without Fermi they would have remained 
undetected for decades or longer. In fact, several 
of them resolved the mystery of what was powering 
unidentified TeV sources. 

For the 4 radio detections we have DM dis- 
tances and can thus convert the fluxes into pseudo- 
luminosities (in units of mjy-kpc 2 ) and these are 
shown in Figure [2] Here we find that one of our pul- 
sars (J2032+4127) is a modestly low luminosity radio 
pulsar while the other three have luminosities of a few 
xlO -2 mjy-kpc 2 , far less luminous than any other 
normal pulsar. This forces us to consider what "radio 
quiet" really means and how pulsars could have such 
a low luminosity. Could we be just scraping the edge 
of the beam or is there some other effect at work? To 
really answer these questions, we will need more ex- 
amples of this low luminosity population and detailed 
modeling of the emission region incorporating geomet- 
ric information from polarization measurements. 

3. Searching for Radio Pulsars in LAT 
Unassociated Sources 

With known pulsars being such a prevalent class of 
GeV 7-ray sources, a complementary strategy is to use 
radio telescopes to search for pulsars powering sources 
discovered in the LAT sky survey. A significant effort 
was put into searches like these during the EGRET 
era, with only modest success [TT]. However, the sen- 
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Figure 2: Pseudo luminosities (5i4ood 2 ) for the four LAT 
blind search pulsars with detected radio pulsations, 
compared to the population of young pulsars in the 
ATNF catalog (blue dots). Three of them are less 
luminous than any previously known young pulsar. 



sitivity and angular resolution of the LAT make these 
searches much more powerful. The sensitivity pro- 
vides many more sources - there were > 600 unasso- 
ciated sources in the First LAT Source Catalog [H H] 
- and the angular resolution is such that most sources 
are localized to a region small enough to be covered by 
a single pointing of a large single-dish radio telescope 
(see Figure [3]). 

The PSC has undertaken a large number of pul- 
sar searches directed at pulsar-like LAT sources. Tar- 
gets are ranked by how much they look like pulsars 
in their spectral and variability characteristics (see 
[5] for the statistical methods used). Initially, tar- 
gets were selected at both high and low Galactic lat- 
itude, but it was quickly realized that the high lat- 
itude searches were highly successful at finding new 
millisecond pulsars, while the low latitude searches 
found only 1 new pulsar associated with a LAT source 
(PSR J2030+3641 [9]) and 3 slow pulsars that are al- 
most certainly chance coincidences unrelated to the 
LAT sources (see Table 

The lack of success at finding young or middle-aged 
pulsars in searches of LAT sources in the plane is dis- 
cussed in [9]. The basic argument is that most de- 
tectable normal radio pulsars in the plane have been 
discovered in large area surveys. Thus, many of the 
pulsar-like 7-ray sources discovered by the LAT near 
the plane may indeed be pulsars, but with the ra- 
dio beams not intersecting our line-of-sight. This is 
supported by the very small fraction (4/35) of blind 
search pulsars from which radio pulsations have been 
discovered. 

The situation at high latitude has been completely 
different. At the end of the EGRET era, there were no 



confirmed 7-ray emitting MSPs. Shortly after launch, 
the LAT demonstrated that MSPs were indeed 7-ray 
emitters, much like their young pulsar brethren [2J. 
Our radio searches of high-latitude LAT sources that 
are non-variable and have curved spectra character- 
istic of pulsars have yielded 43 new MSPs (42 were 
discovered by the PSC and one by a non-PSC group 
also targeting LAT source; see Table [TT] and r igure Ion . 

Also, in contrast to the situation for young pulsars 
in the plane, the sky at high latitudes has mostly 
not been surveyed with good sensitivity to millisec- 
ond pulsars. This is changing with several new radio 
surveys designed to detect MSPs (e.g. HTRU, GBT 
Drift, GBNCC). In several cases, new MSPs discov- 
ered in these surveys have turned out to be LAT pul- 
sars as well |19) . and more are sure to come. 

As the newly-discovered pulsars are timed over a 
year or so after discovery, we obtain orbital parame- 
ters and precise positions enabling us to fold the LAT 
7-ray data to look for pulsations that confirm the asso- 
ciation with the 7-ray source. In nearly all cases where 
phase-coherent timing solutions have been obtained 
the associations have been confirmed (as marked in 
Table |II|), but in one case (PSR J1103-5403) the pul- 
sar was not associated with the LAT source but just 
in the radio beam by chance. This is likely the case for 
more of the new discoveries with the GMRT, owing to 
its very large beam size. 

One striking feature of this new population of MSPs 
discovered in 7-ray directed searches is the fraction 
in 'black widow' (pulsars in orbit around very low 
mass companions and usually eclipsing) and 'redback' 
(eclipsing binaries with a main sequence companion; 
see [25]) systems. At the time of the Fermi launch, 
there were about 60 MSPs known outside of the glob- 
ular cluster system. These included 3 black widows 
and 1 redback. Among the 43 new pulsars found in 
these searches, there are at least 10 black widows and 
4 redbacks (additions to these numbers are possible 
because some of the pulsars have not yet had their 
orbital parameters determined). One possible expla- 
nation for this is that it is a selection effect, where the 
energetic MSPs powering 7-ray sources are more likely 
to be in these types of interacting binary systems, or 
possibly that some of the 7-rays are coming from the 
interaction region itself (intra-binary shock). 

Another indication that this population is some- 
what different than the population discovered in large 
area radio surveys is the period distribution (see Fig- 
ure [4| , which is much more strongly clustered at the 
short period end of the distribution than the general 
field MSP population (defined here as the 63 MSPs 
(P < 16 ms) in the ATNF catalog with discovery 
dates before 2010 and not associated with globular 
clusters) . The median periods of the two samples are 
3.16 ms for the PSC MSPs and 4.62 ms for the previ- 
ously known MSP population, and the K-S test give 
a probability of 1.3 x 10 -4 that they are drawn from 
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the same underlying distribution. All of this supports 
earlier studies that find that the spin periods in eclips- 
ing systems are preferentially faster than non-eclipsing 
systems [14"] . 

This large number of new MSPs also provides 
many potential new additions to pulsar timing array 
projects. Most of the eclipsing systems exhibit or- 
bital period evolution and other effects that keep them 
from being among the best MSP clocks, but seven of 
the new systems are being timed by the NANOgrav 1 
project to evaluate their potential contribution to the 
detection of gravitational waves. 

So, while these new MSPs are all radio pulsars that 
could eventually have been discovered in undirected 
radio surveys, Fermi has enabled us to discover them 
in a much more efficient manner and much earlier than 
they would have been, allowing for detailed follow-up 
observations and potentially speeding up the quest for 
the direct detection of gravitational waves. With dif- 
ferent biases than the radio surveys, we are helping to 
get a more complete picture of the Galactic MSP pop- 
ulation and finding many systems that will improve 
our understanding of MSP formation and evolution. 
One lesson has been made very clear in our searches. 
Each good candidate LAT source must be searched 
multiple times before declaring it to be devoid of a 
detectable radio pulsar. In many cases the eclipses or 
strong interstellar scintillation prevent detection in a 
substantial fraction of the observations of a particular 
source. Recently we have been focusing on making 
sure that the best candidate sources are observed 3, 
4, or more times to account for this. 

But, what about the pulsar- like high latitude 
sources where no radio pulsations are ever discov- 
ered? These could be a population of radio quiet (as 
seen from Earth) MSPs. Direct searches for millisec- 
ond pulsations in the LAT data are extremely com- 
putationally challenging, even for the ~20% of MSPs 
that are isolated PUdl]. For tlie binary MSPs, these 
searches are not currently tractable. However, given 
the large fraction of black widow and redback systems 
among the LAT-detected MSPs, another approach is 
promising. Optical and X-ray studies of the fields of 
the bright pulsar-like LAT source OFGL J2339.8-0530 
revealed a likely counterpart that is almost certainly 
an eclipsing binary MSP system [HJ [30]- With the 
orbital parameters and precise position determined 
from optical observations, a search for 7-ray pulsa- 
tions from the pulsar in this system may be feasible. 
X-ray and optical studies of the error boxes of other 
LAT sources may reveal more examples of this kind 
of system. 



Table I Young Pulsars Discovered in PSC Searches 



PSR 


Period 


DM 


LAT PSR? Reference 




(ms) 


(pc cm -3 ] 




J1604-44 


1389.20 


176 


No [18] 


J2030+3641 


200.13 


247 


Yes [5] 


J1400-56 


410.70 


123 


No [20] 


J1203-62 


393.09 


285 


No [20] 




1FGL Localization 



30 

conf) 



1 http : //nanogr av . org 



Figure 3: Cumulative number of sources by 95% 
confidence localization in the 1FGL catalog. The FWHM 
beam widths of several telescope/frequency combinations 
are labeled. Notice that Parkes at 1400 MHz and GBT 
at 820 MHz can cover 800-1000 of the 1400 sources in a 
single pointing. By going to 350 MHz, even the 100-m 
GBT can cover the full error circle for nearly every 
source in the 1FGL catalog. The GMRT (the most recent 
addition to the PSC) has an extremely large beam width 
of 40' at 610 MHz and 80' at 325 MHz! 
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Table II Millisecond pulsars discovered in radio pulsation searches of LAT unassociated sources 



PSR 


Telescope 


LAT Source 




P 


DM 


D 




-^^c,min 


Notes 


Refs 










H 


(pc cm ' 


! ) 


(kpc) 


(d) 


(M @ ) 






J0023+0923 


GBT-350 


1FGLJ0023. 


.5+0930 


3. 


.05 


14 


.3 





.7 


0.14 


0.017 


BW, 7 


m 


J0101-6422 


Parkes 


1FGLJ0101. 


.0-6423 


2 


.57 


12 


.0 





.6 


1.78 


0.16 


7 




J0102+4839 


GBT-350 


1FGLJ0103. 


.1+4840 


2 


.96 


53 


.5 


2 


.3 


1.67 


0.18 


7 


m 


J0307+7443 


GBT-350 


1FGLJ0308. 


.6+7442 


3 


.16 


6 


.4 





.6 


36.98 


0.24 


7 


.16 


J0340+4130 


GBT-350 


1FGLJ0340. 


.4+4130 


3. 


.30 


49 


.6 


1 


.8 


Isolated 




7 


nu 


J0533+67 


GBT-820 


2FGLJ0533. 


.9+6759 


4 


.39 


57 


.4 


2 


.4 


Isolated 








J0605+37 


GBT-820 


2FGLJ0605. 


.3+3758 


2. 


.73 


21 


.0 





.7 


55.6 


0.18 






J0614-3329 


GBT-820 


1FGLJ0614. 


.1-3328 


3. 


.15 


37 


.0 


1 


.9 


53.6 


0.28 


7 


m 


J0621+25 


GBT-820 


2FGLJ0621. 


.2+2508 


2 


.72 


83 


.6 


2 


.3 


TBD 








J1103-5403 


Parkes 


1FGLJ1103. 


.9-5355 


3. 


.39 


103 


.9 


2 


.5 


Isolated 




NA 


m 


J1124-3653 


GBT-350 


1FGLJ1124 


.4-3654 


2 


.41 


44 


.9 


1 


.7 


0.23 


0.027 


BW, 7 




J1124-36 


GMRT-325 


1FGLJ1124. 


.4-3654 


5. 


.55 


45 


.1 


1 


.7 


TBD 




NA 




J1142+0119 


GBT-820 


1FGLJ1142. 


.7+0127 


5. 


.07 


19 


.2 





.9 


1.58 


0.15 






J1207-5050 


GMRT-610 


1FGLJ1207. 


.0-5055 


4 


.84 


50 


.6 


1 


.5 


Isolated? 








J1231-1411 


GBT-820 


1FGLJ1231. 


.1-1410 


3 


.68 


8. 


.1 





.4 


1.86 


0.19 


7 


mi 


J1301+0833 


GBT-820 


1FGLJ1301. 


.8+0837 


1 


.84 


13 


.2 





.7 


0.27 


0.024 


BW 




J1302-32 


GBT-350 


1FGLJ1302. 


.3-3255 


3. 


.77 


26 


.2 


1 


.0 


0.78 


0.15 




nu 


J1312+0051 


GBT-820 


1FGLJ1312. 


.6+0048 


4 


.23 


15 


.3 





.8 


38.5 


0.18 


7 




J1514-4946 


Parkes 


1FGLJ1514. 


.1-4945 


3 


.59 


30 


.9 





.9 


1.92 


0.17 


7 


m 


J1536-4948 


GMRT-325 


1FGLJ1536. 


.5-4949 


3. 


.08 


38 


.0 


1 


.8 


TBD 








J1544+4937 


GMRT-610 


18M3037 




2 


.16 


23 


.2 


1 


.2 


0.117 


0.018 


BW 


m 


J1551-06§ 


GBT-350 


1FGLJ1549. 


.7-0659 


7 


.09 


21 


.6 


1 


.0 


5.21 


0.20 




m 


J1628-3205 


GBT-820 


1FGLJ1627 


.8-3204 


3. 


.21 


42 


.1 


1 


.2 


0.21 


0.16 


RB 




J1630+37 


GBT-820 


2FGLJ1630. 


.3+3732 


3. 


.32 


14 


.1 





.9 


12.5 


0.16 






J1646-2142 


GMRT-325 


1FGLJ1645 


0-2155c 


5 


.85 


29 


.8 


1 


.1 


23 


TBD 






J1658-5324 


Parkes 


1FGLJ1658 


.8-5317 


2 


.43 


30 


.8 





.9 


Isolated 




7 


m 


J1745+1017 


Effelsberg 


1FGLJ1745. 


.5+1018 


2 


.65 


23 


.9 


1 


.3 


0.73 


0.014 


BW, 7 


m 


J1747-4036 


Parkes 


1FGLJ1747. 


.4-4036 


1. 


.64 


153 


.0 


3 


.4 


Isolated 




7 




J1810+1744 


GBT-350 


1FGLJ1810. 


.3+1741 


1. 


.66 


39 


.6 


2 


.0 


0.15 


0.045 


BW, 7 


m 


J1816+4510 


GBT-350 


1FGLJ1816. 


.7+4509 


3. 


.19 


38 


.9 


2 


.4 


0.36 


0.16 


t, RB- 7 


m 


J1828+0625 


GMRT-325 


1FGLJ1830. 


.1+0618 


3. 


.63 


22 


.4 


1 


.2 


6.0 


TBD 






J1858-2216 


GBT-820 


1FGLJ1858. 


.1-2218 


2 


.38 


26 


.6 





.9 


46.1 


0.22 


7 




J1902-5105 


Parkes 


1FGLJ1902. 


.0-5110 


1 


.74 


36 


.3 


1 


.2 


2.01 


0.16 


7 


m 


J1902-70 


Parkes 


2FGLJ1902. 


.7-7053 


3. 


.60 


19 


.5 





.8 


TBD 








J2017+0603 


Nangay 


1FGLJ2017. 


.3+0603 


2 


.90 


23 


.9 


1 


.6 


2.2 


0.18 


7 


nni 


J2043+1711 


Nangay 


1FGLJ2043 


.2+1709 


2 


.38 


20 


.7 


1 


.8 


1.48 


0.17 


7 


m 


J2047+1053 


GBT-820 


1FGLJ2047. 


.6+1055 


4 


.29 


34 


.6 


2 


.0 


0.12 


0.036 


BW, 7 




J2129-0429 


GBT-350 


1FGLJ2129. 


.8-0427 


7 


.62 


16 


.9 





.9 


0.64 


0.37 


RB 


ED 


J2214+3000 


GBT-820 


1FGLJ2214. 


.8+3002 


3. 


.12 


22 


.5 


1 


.5 


0.42 


0.014 


BW, 7 


26, 


J2215+5135 


GBT-350 


1FGLJ2216 


.1+5139 


2 


.61 


69 


.2 


3 


.0 


0.17 


0.21 


RB, 7 


nu 


J2234+0944 


Parkes 


1FGLJ2234. 


.8+0944 


3 


.63 


17 


.8 


1 


.0 


0.42 


0.015 


BW, 7 


m 


J2241-5236 


Parkes 


1FGLJ2241. 


.9-5236 


2 


.19 


11 


.5 





.5 


0.14 


0.012 


BW, 7 


m 


J2302+4442 


Nangay 


1FGLJ2302. 


.8+4443 


5 


.19 


13 


.8 


1 


.2 


125.9 


0.3 


7 


m 



Notes: The number after the telescope name indicates the observing frequency in MHz; no number means 1400 MHz. 
§= This pulsar was previously published as J1549— 06 |16l 123] before the precise position was determined. 
BW = black widow, RB = redback, NA = known to not be associated with the target LAT source, 
f = Discovered by GBNCC team in analysis of beams coincident with LAT sources ,34 , 7 = 7-ray pulsations detected 
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Figure 4: Histograms of pulse periods of MSPs 
discovered before Fermi and those discovered in PSC 
searches of LAT unassociated sources. The selection for 
7-ray emission selects for faster MSPs than large area 
radio surveys. 
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Figure 5: Fermi LAT sky map in Galactic coordinates showing the locations of the 43 MSPs discovered in searches of 
LAT unassociated source error circles. The symbols are color coded by the discovery telescope: blue for GBT, red for 
Parkes, green for GMRT, magenta for Nangay, black for Effelsberg. One GMRT discovery (J1124-36) is not shown since 
it overlaps in position with J1124-3653. The one square red symbol for J1103-5403 indicates that it was a serendipitous 
discovery of an MSP not associated with the LAT source [TH]. One of the GBT pulsars (J1816+4510) was discovered 
by the GBNCC team (PI: S. Ransom) in an analysis of survey beams coincident with LAT catalog sources [34] . 
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